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Abstract: Local variety of lime ( Lemun tsami) was subjected to 
quasi-static compressive loading using Instron Universal 
Testing Machine at three levels of time after harvest; freshly 
harvested , one and two weeks after harvest Creep parameters 
were determined by fitting the experimental data to the 
Burger’s Model and the obtained model equation for freshly 
harvested was: J(t ) = 1. 843 + 2. 59[ 1 - e ( ' t/2 - m> ]+ — , 

R 2 = 0.98 . The instantaneous creep compliance (J 0 ) was; 

1.843, 3.137 and 4.116 1/MPa for freshly harvested, one and 
two weeks after harvest respectively while the corresponding 
instantaneous elasticity (E 0 ) was 0.543, 0.319 and 0.243 MPa 
respectively. Compliance increased with increase in time after 
harvest while elasticity decreased with time after harvest. The 
obtained free viscosity (q 0 ) was 1160.093, 808.407 and 808.859 
MPas for freshly harvested, one and two weeks after harvest 
respectively. This was found to decrease with increase in time 
after harvest. 

Keywords: Lime (Lemun tsami), Creep, compliance function, 
Retardation time, Instantaneous Viscosity 

I. Introduction 

Lime (Citrus aurantifolia) is a citrus species with a globose fruit. 
It is the smallest as compared to other citrus species with a 
thinner rind. It has many uses; the juice is used to flavor food 
and cola drinks, preparation of pickles and sauces. Although the 
oil from its rind can irritate the skin because of its high acidic 
content, it has a lot of medicinal uses; it is antiseptic, antiviral, 
bactericidal and disinfectant ( Wikipedia, the free encyclopedia) 
Despite the importance of lime, its production is said to be on 
the decline worldwide due to incidence of citrus canker 
(Hayley, 2014; Thomas, 2010). In Nigeria, there is little 
information on the engineering properties of the cultivar, this 
may be due to the fact that lime is used on a small scale as 
compared to other citrus cultivars as affirmed by Aiyelaagbe, et 
al f . (1996). 

Lime is susceptible to injuries such as bruises, cuts, and puncture 
during harvesting, handling, packaging, transportation and 
storage. The situation is made worst by the transportation of this 
fruits in bulk in wagons under harsh environmental conditions 



and on bad roads. According to Nabil (2013), the knowledge of 
viscoelastic properties of fruits and vegetables do not only 
minimizing damage during harvesting, transportation and 
storage but also help in understanding physiological changes that 
take place in the fruit during growth, maturation, ripening and 
storage after harvest as well as assessment of product quality. 
Also, the knowledge of viscoelastic property of fruits is 
important in quality evaluation, equipment development and 
product processing. 

Viscoelastic properties of fruits and vegetables are characterized 
by two distinctive phenomena; creep and stress relaxation. The 
major creep parameters include: Creep compliance function J(t), 
a measure of relative ease with which a material can be 
deformed under a static load, retardation time (T ret ), the time 
taken for the material to deform to ( 1 — - ) or about 63.2 % of 
the total deformation at rupture point under static load, and free 
viscosity (q). 

The objectives of this work are to determine: (i) Creep 
compliance function J(t), retardation time (T ret ), and Viscosity 
(r|), at three levels of time after harvest (freshly harvested, 1 
week and 2 weeks after harvest) respectively; and (ii) effects of 
time after harvest on these parameters. 

Mathematical Model 

The rheological model to represent the creep behavior of food 
and biological materials are described by the four element 
Burger’s model (Mohsenin, 1986 and Kelly, 2011). In fig. 1, 
when a static load is suddenly imposed on the material, spring E 0 
is stretched suddenly by an amount (O - A) which is e G = o 0 /E 0 
After this initial strain, e 0? creep starts at a high rate but gradually 
slows down due to dashpot r| 0 resulting in retarded elastic 
behavior (A-B). After time, t, (at the end of loading), when the 
imposed load is removed spring E 0 snaps back to its original 
state but spring E r cannot contract to its original state 
instantaneously because of rji and rj 0 but does so gradually. 

When the load is removed, the spring E r slowly forces the 
plunger of dashpot r|i back to its original position; however, 
since no force is acting on the dashpot r| 0 upon load removal, this 
element retains a non-recoverable displacement representing a 
permanent deformation in the material ( Mohsenin, 1986). 
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Fig.l: Mechanical Analogy of Four Element Burger’s Model for 
a Typical Biological Material 



The mathematical expression for Burger’s Model is given by 
equation (1) ( Mohsenin, 1986; Abd El-Maksoud, et al . , 2009; 
Kelly, 2011; Meyers, 2011) 



e(t) 



° + JL +^ [1 _ e (-^«) ]+ ^ 

E 0 Er rjo r| 0 



( 1 ) 



II. Materials and Methodology 

The local cultivar of lime ( Lemun tsami) was obtained from 
Kaura Citrus Farm in Toto Local Government Area of Nasarawa 
State, North Central Nigeria. Four trees in plots of trees typical 
of the variety were selected from which fruits were harvested for 
the tests. 



Some fruits were carefully handpicked from the trees while 
others were chipped off the three with a knife leaving a stalk 10 
- 12 cm long and leaves removed (Coppock, et al, 1969 ); this is 
to maintain some level of physiological freshness for tests 
concerning freshly harvested. The fruits were kept cooled in a 
fruit shed by water spray while harvesting was going on; at the 
end of harvest (between 1.00 - 2.00 pm), they were then packed 
in cardboard boxes. The bottoms of these boxes were lined with 
foam to minimize mechanical injuries and sides perforated to 
reduce temperature and ethylene build up (Tabatabaekoloor, 
2012). 

They were transported the same day to Advanced Materials 
Laboratory of the Engineering Materials Development Institute 
(EMDI), KM 4, Ondo Road, Akure, Ondo State, Southwest 
Nigeria and stored in a cool room maintained at about 5 °C, and 
87 % relative humidity immediately upon arrival at about 8.30 
pm. Tests for freshly harvested was conducted at 7.30 am the 
following day (about 11 hours after harvest) while other tests 
were conducted after 7 and 14 days respectively. 



mean dimensions; major, intermediate and minor diameter were 
found to be 3.876 ± 0.409, 3.551 ± 0.291 and 3.512 ± cm 
respectively while geometric mean diameter and sphericity were 
respectively 3.640 + 0.304 cm and 0.942 ± 0.042. 



Creep Tests 

Preliminary tests were first conducted on the samples to estimate 
the mean constant force that would produce deformation of 
about 63.2 % (equivalent of 1 — of the total deformation at 
rupture point for each of the treatments (Mohsenin, 1986). 

This was done by grouping the fruits into two, each group 
having nearly the same physical characteristics (size and 
sphericity). Surface moisture on the specimen was cleaned, and 
then placed axially in the Instron Universal Testing Machine 
(Model 3369, No. K334; 50 kN capacity) under parallel steel flat 
plate (Plate 1) and loaded to rupture point. Each test was 
replicated five times for each group at the rate of 1 mm/s and the 
mean force obtained. The same procedure was followed for all 
the treatments. 



Having obtained the mean force and the corresponding 
deformation at rupture point for all the treatments, 63.2 % of the 
constant force obtained in the preliminary test was imputed in to 
the machine as static load for each group. Each specimen was 
loaded at 1 mm/s ( with the 63.2 % mean load obtained in the 
preliminary test) for 920 s, the strain with time was 
automatically captured by the micro -computer of the machine ( 
during loading only) and the corresponding constant stress (o 0 ) 
of each specimen calculated using equation (2) according to 
Jatuphong, et al., (2008) 



Constant force (. F 0 ) A T 

o 0 = •• •• (2) Although the 

Final Contact Area (A 0 ) "" 

area of contact of test specimen varies during loading, the 
contact area registered by the machine at the end of loading was 
used (Jatuphong, et al., 2008). 



The test was replicated five times for each group for all the 
treatments and the mean strain at each time interval was 
obtained. 



Creep Compliance. 

Creep compliance ( J ) is the ratio of strain to stress ( - ) or the 

o 

reciprocal of modulus (-), thus diving equation (1) by the 
calculated constant stress (o 0 ) yields equation (3) which is the 
creep compliance function J(t). 

J(f)=Jo+Jr[ l-e ( -^ re,) ]+^ .. (3) 



Size and Shape 

Dimensions of each specimen were determined on the three 
mutually perpendicular axes using a digital vernier caliper 
reading to 0.01mm. The fruits were characterized in terms of 



The mean creep compliance of each treatment is presented in 
Table 1. 

III. Results and Discussions 
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Table 1: 


Mean Creep Compliance (1/Mpa) 


Time (s) 


Time after Harvest (Days) 




0 


7 


14 


0 


1.843 


2.963 


3.939 


5 


1.980 


3.137 


4.116 


10 


2.106 


3.302 


4.284 


15 


2.222 


3.450 


4.432 


20 


2.317 


3.579 


4.565 


25 


2.399 


3.697 


4.698 


30 


2.478 


3.807 


4.802 


35 


2.556 


3.909 


4.902 


40 


2.614 


4.007 


5.002 


45 


2.680 


4.076 


5.082 


50 


2.722 


4.166 


5.158 


55 


2.770 


4.208 


5.218 


60 


2.815 


4.274 


5.283 


70 


2.854 


4.324 


5.329 


80 


2.888 


4.372 


5.410 


90 


2.921 


4.423 


5.419 


100 


2.939 


4.456 


5.474 


110 


2.975 


4.492 


5.490 


120 


2.998 


4.517 


5.514 


160 


3.009 


4.553 


5.540 


200 


3.037 


4.581 


5.568 


240 


3.055 


4.608 


5.592 


280 


3.072 


4.630 


5.600 


320 


3.082 


4.652 


5.618 


370 


3.096 


4.668 


5.631 


420 


3.102 


4.684 


5.658 


470 


3.108 


4.701 


5.673 


520 


3.116 


4.714 


5.693 


570 


3.122 


4.734 


5.700 


620 


3.132 


4.736 


5.712 


670 


3.138 


4.755 


5.726 


720 


3.140 


4.760 


5.733 


770 


3.142 


4.760 


5.793 


820 


3.144 


4.769 


5.796 


870 


3.145 


4.770 


5.797 


920 


3.147 


4.772 


5.800 



Creep Compliance - Time Curves 
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Fig. 2: Creep Compliance of Freshly Harvested Time 
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Fig. 3: Variation of Creep compliance with time after harvest 

Figs. 2 is the creep compliance - time curve of freshly harvested 
Fime while Fig. 3 is the variation of creep compliance with time 
after harvest. 

Analysis of Creep Parameters from Creep Compliance - 
Time Curve 

The Model parameters viz., instantaneous Compliance (, J 0 ), 
instantaneous elasticity ( E 0 ), Viscosity of the free dashpot ( f ] Q ), 
retarded compliance (Jr), retarded elasticity (Er), retarded 
viscosity (i/ r ), and retardation time (T ret ) were obtained from the 
Compliance data as illustrated by Fig. 4 following steps (i) - 
(vii) 




Fig.4: Graphical illustration for calculation of creep parameters 

i) J 0 = Instantaneous response (O - A). For freshly 
harvested, J 0 = 1.843 MPa' 1 
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ii) From J 0 = — , E 0 = 0.543 MPa 

E 0 

iii) The slope of the straight line portion of the curve (B - 
C) equals — while the intercept gives Jr. This was 

obtained by linear regression analysis using IBM ® 
SPSS ® Statistics, Version 20.0 
For freshly harvested, the slope is 0.000862, therefore tj a = 
1160.093MPas, the intercept is 2.592, therefore Jr = 
2.592 MPa 1 

iv) From J, = — , E, = 0.386 MPa 

Er 

v) However, retardation time ( T ret ) was determined from 
the exponential portion of the curve (A - B) (Gorji, et 
al., 2010) whose equation is of the form: 

J(t) — Jo = Jr | 




Fig. 6: Variation of Instantaneous and Retarded Compliance 
With Time After Harvest 



Re-arranging and taking log on both sides gives 
equation (4) Ln [l - ^— ] = ^ ( 4 ) 

Inserting appropriate values, T ret = 2.933 sec. 

vi) From T ret = r| r / E r , r| r = 0.1 13 MPas. 

vii) Putting these model parameters in equation (3) gives 
equation (5) for freshly harvested cultivar 

J (t) = 1.843 + 2.59[ 1 - e M/2 ' 933) ]+ — - — , R 2 = 0.98 .. 

L J 1160.093’ 

.. (5) 

This procedure was followed and the parameters obtained 
for all the treatments are presented in Table 2. 



Table 2: Creep parameters for Time 


Parameters 


Time after Harvest 




Fresh 


1 week 


2 weeks 


Jo (1/MPa) 


1.843 


3.137 


4.116 


Eo(MPa) 


0.543 


0.319 


0.243 


Jr (1/MPa) 


2.592 


3.965 


4.956 


Er(MPa) 


0.386 


0.252 


0.202 


rjr (MPas) 


0.113 


0.189 


0.886 


Tret (s) 


2.933 


3.75 


4.386 


rj 0 (Mpas) 


1160.093 


808.407 


808.859 
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Fig. 7: Variation of Viscosity (q 0 ) with Time after Harvest 
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Fig. 8: Variation of Retardation Time (T ret ) With Time after 
Harvest 




Fig. 5: Variation of Instantaneous and Retarded Elasticity with 

time after Harvest Fig. 9: Variation of Creep Parameters with time after harvest. 
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Fig.2 shows that creep compliance has an initial exponential 
increase with time but levels off as time increases to infinity. 
This means that as the static load is suddenly applied, it is 
initially borne by both the elastic and viscous components to a 
point where the elastic component ceases; where the viscous 
component dominates resulting in the curve leveling off as 
affirmed by Gorji et al, (2010) for apple and Mohsenin (1986) 
for most fruits and vegetables. Fig. 3 shows the variation of 
compliance curve with time after harvest. 

Instantaneous elasticity was found to be 0.543, 0.319 and 0.243 
MPa for freshly harvested, 7 and 14 days after harvest 
respectively; while retarded elasticity was 0.386 0.252 and 0.202 
MPa for the corresponding time after harvest. However, Gorji et 
al, (2010) obtained from prepared samples of two varieties of 
Apple ( Golab and Shafi Abadi) loaded at 25.4 mm/min an 
instantaneous elasticity of 714.28 and 400 MPa respectively, 
while the retarded elasticity was 7.7 and 7 MPa repectively. 

The differences between the values of instantaneous elasticity 
obtained for Fime as compared to Apple may be due to the 
differences in the internal structures of the fruits. According to 
Panmanas and Charoonpong (2012), unlike most fruits whose 
internal structures are homogeneous, citrus fruits have high 
compartmentalized internal structure; juice sacs (which contain 
juice) are enclosed in a segment covered by a tough lamella, 
these segments are enclosed in a rind made up of Albedo (which 
has white spongy texture) and Flavedo (greenish part dotted with 
oil glands. These structures give the fruits its high elastic 
property. 

In fig. 5, both elasticity (E 0 , and E r ,) decreased with time after 
harvest. This is similar with what Nabil (2013) and Ayman, et 
al, (2012) observed in Tomato fruits and Pears respectively 
during storage; that both instantaneous and retarded elasticity of 
these fruits decreased with time of storage. Nabil (2013) stated 
that the obtained values of both instantaneous and retarded 
elasticity for Tomato fruits were found to be inversely 
proportional with storage time. This phenomenon may probably 
due to ripening process. According to Naoki and Donald (1997), 
elasticity in fruit tissues change as a function of fruit ripening; 
during ripening, the polysaccharide constituents (cellulose, 
pectin) in cell walls are degraded by wall hydrolyzing enzymes 
thus softening the textural properties of the fruits. 

Instantaneous Creep Compliance (, J 0 ) was found to be 1.843, 
3.137, and 4.116 1/MPa for freshly harvested, 7 and 14 days 
after harvest respectively; while retarded Compliance (, J r ) was 
2.592, 3.965 and 4.956 1/MPa for the corresponding time after 
harvest. 

I n Fig. 6, the following were observed: (i) the values of 
instantaneous compliance (, J 0 ), was greater than that of retarded 
compliance (, J r ) for the period under review. This is because the 
J 0 represents the elastic component (with higher relative ease to 
deform and recover) while the J r is a combination of viscous and 
elastic component. (ii) Both instantaneous and retarded 
compliance increased with time after harvest. This may probably 
be due to ripening process; according to Naoki and Donald 



(1997), the polysaccharide constituents (cellulose, pectin) in cell 
walls are degraded by wall hydrolyzing enzymes during ripening 
thus softening the textural properties of the fruits. 

Instantaneous viscosity (rj a ) was found to be 1160.093, 808.407 
and 808.859 MPas for Freshly harvested, 7 and 14 days after 
harvest respectively. However, for two varieties of Apple ( Golab 
and Shafi Abadi ) loaded at 25.4 mm/min; Gorji et al, (2010) 
obtained rj 0 = 2000 MPas for both varieties when freshly 
harvested while for Tomato fruits, Nabil (2013) obtained values 
of 1500 and 2500 N.min/mm for axial and longitudinal loading 
positions respectively. These differences in values of the three 
crops may not be unconnected with differences in their internal 
structures as earlier explained. In fig. 7, the instantaneous 
Viscosity ( t ] 0 ) decreased with increase in time after harvest, this 
is in accordance with what Nabil (2013) obtained for Tomato 
fruits. 

The retardation time were found to be 2.933, 3.750, 4.386 s for 
freshly harvested, 7 and 14 days after harvest respectively. 
However, Gorji et al, (2010) obtained retardation values of 12 
and 15 s for varieties of apple, Golab and Shafi Abadi 
respectively. This implies that lemun tsami has the ability to 
deform within a short time period and also recovers at the same 
time period upon removal of the imposed load. In comparison, 
lemun tsami is more elastic than Apple. In fig. 8, retardation 
time ( T ret ) increased as time after harvest increases. Fig. 9 
shows a graphical variation of creep parameters with time after 
harvest. 

IV. Conclusion 

From the results obtained, Lemun tsami is best handled when 
fresh; for at this state, it has relatively less tendency to deform 
under dead load (Jo = 1.843 1/MPa), more elastic ( E 0 = 0.543 
MPa) and can dissipate internally imposed strain at a relatively 
short time ( T ret = 2. 933 s). 
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Fig. 1 : Lime undergoing quasi-static compression 
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